Effects of molecular adsorption on optical losses of Ag (111) surface. 



A. V. GavrilenkoQ C. S. McKinney, and V. I. Gavrilenko 
Center for Materials Research, Norfolk State University, 700 Park Ave, Norfolk, VA 23504 

The first principles density functional theory (DFT) is applied to study effects of molecular ad- 
sorption on optical losses of silver (111) surface. The ground states of the systems including water, 
methanol, and ethanol molecules adsorbed on Ag (111) surface were obtained by the total energy 
minimization method within the local density approximation (LDA). Optical functions were calcu- 
lated within the Random Phase Approximation (RPA) approach. Contribution of the surface states 
to optical losses was studied by calculations of the dielectric function of bare Ag (111) surface. 
Substantial modifications of the real and imaginary parts of the dielectric functions spectra in the 
near infrared and visible spectral regions, caused by surface states and molecular adsorption, were 
obtained. The results are discussed in comparison with available experimental data. 



I. INTRODUCTION 

The effects of metallic surface states and molecular ad- 
sorption on the atomic surface structure and electronic 
properties of transition metals is extensively discussed in 
literature (see e.g. Refs.[TMI and references therein). The 
silver based metallic nanostructures are actively stud- 
ied for numerous applications of modern nanotechnol- 
ogy such as nanophotonics and nanoplasmoniaffl It has 
been demonstrated before that electronic surface states 
on metallic surfaces substantially contribute to the bond- 
ing of physisorbed ato ms an d molecules^ and to their op- 
tical absorption spectra^'. The semi-classical approach 
developed by PerssorP using Drude model and modified 
by Pinchuk et alP clearly demonstrated increasingly im- 
portant contributions of the surfaces states on optical ab- 
sorption spectra of metallic nanoparticles with descrease 
of their size. 

Reduction of optical losses in nanostructured materi- 
als is one of the most important issues affecting success- 
ful engineering of different optoelectronics photonics de- 
vices. Optical losses of materials are determined through 
the imaginary part of the dielectric function (Im(e))' 7 . 
Though losses could be compensated through external 
gain, another possibility is the modification of electronic 
structure of metallic nanoparticles by chemical adsorp- 
tion on the surfaces. However, this approach requires 
detailed understanding of electronic processes accompa- 
nied by chemical adsorption on a microscopic scale. 

Recently using the first principle modeling we pre- 
dicted the effect of the surface states on the dielec- 
tric function spectra (s(oj)) of the nanometer-thick silver 
filmiP. It has been shown that a decrease of the silver- 
slab thicknesses resulted in substantial increase of the 
Im(e) due to the surface state contributions. Indepen- 
dently, these predictions have been confirmed experimen- 
tally by Drachev et alP by simultaneous measurements 
of reflectance and transmittance of differently shaped sil- 
ver nanostructures. They show that the dramatic in- 
crease of Im[e(uj)] for both TE and TM light polariza- 
tions is due to the decrease of characteristic dimensions 
of the silver-based nanostructures; agreeing well with our 
predictions^. 



In the present work, we study physical and chemical 
processes affecting optical losses of silver surfaces that 
substantially contribute to the optics of metamaterials 
(the artificially fabricated nano- materials^). The adsorp- 
tion of water, ethanol, and methanol on silver nano- 
films is studied by first principle density functional theory 
(DFT). These systems are w idely used in modern nan- 
otechnology and nano-physics^^. Contribution of the 
surface states to optical losses is studied by calculations 
of the imaginary part of the dielectric function for metal- 
lic slabs of different thickness. We obtained substantial 
modifications of the optical function of metallic nano- 
slabs in the near infrared and visible spectral regions that 
are discussed in comparison with literature. 

II. METHOD 

Typical metallic nanoparticles are much bigger in size 
than the molecules studied. Consequently the interac- 
tions of small molecules with the nano-particles can be 
studied using widely used model of the surface physics: 
molecule on an atomic slab^n Equilibrium atomic struc- 
tures of silver surfaces with organic molecules are ob- 
tained from the total energy minimization method within 
DFT using ab initio (gradient-corrected) functional d 11 * 12 !. 
The basis functions are given numerically as values on 
an atomic-centered spherical-polar mesh, rather than as 
analytical functions (i.e. Gaussian orbitals), as imple- 
mented in the DMol 3 package 1 -. We employ the super- 
cell approach to model the surface; our unit cell is a slab 
consisting of seven monolayers of Ag (111) — V3 x \/3, 
of which the top three layers were allowed to relax. A 
Monkhorst Pack mesh of 6 x 6 x 1 k-points was used to 
sample the Brillouin zone. Explicit tests were performed 
to verify that the results were converged for the param- 
eters given above. 

The electronic structure and optical properties of 
the system were then calculated using ab initio norm- 
conserving (NC) pseudopotentials (PP), as implemented 
in the CASTEP package^. The eigen value problem was 
solved using the first principle PP method with NC-PF'S'. 
Convergence of the results has been carefully checked by 
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choosing energy cut-off values up to lOOOeV. The cut- 
off energy of 550eV was chosen for NC-PP to ensure 
reliable convergence of predicted optical data. Choice 
of pseudopotentials (ultrasoft or norm-conserving) does 
not significantly affect the results. Exchange and cor- 
relation (XC) interaction is modeled using the Perdew- 
Burke-Ernzerhof formalism^. A Monkhorst Pack grid of 
6x6x1 was used to sample the Brillouin zone, just like 
for DMol 3 , yielding well converged results. 

Calculated self-consistent eigen energies and eigen 
functions are used as inputs for optical calculations. 
Within its penetration depth beneath the surface, the 
incident light field (E(u>)), at frequency u), induces a lin- 
ear optical susceptibility function x ■ 

P i {w) = X ^E j {u) (1) 

The function x is then used to determine the com- 
plex dielectric function (e = Ei + 1E2) through the linear 
optical susceptibility: e = 1 + 4ttx^ ■ Within the Ran- 
dom Phase Approximation (RPA) the function e is given 
bjEl 



III. RESULTS AND DISCUSSION 

Realistic modeling of optical functions of solid surfaces 
and interfaces still remains challenging for the first prin- 
ciple theories^. As a first step it requires reliable re- 
production of the ground state of the material (equilib- 
rium atomic geometry) by minimizing the free energy. 
As a next step optical properties could be calculated us- 
ing the optimized geometry as an input. Description of 
excited states, which could be in good agreement with 
experiment, normally requires inclusion of local field, 
many-body (excit ons) e ffects, and probably other non- 
local contribution d 18 * 1 ^ . This makes theory much more 
complicated than frequently used independent particles 
approach (or Random Phase Approximation, RPA). In 
this work we focus on relative changes of electron en- 
ergy structure and optical response caused by molecular 
adsorption on Ag (111) surface, which could be realis- 
tically d escrib ed without inclusion of many-body effects 
in optic3^22l. As such we chose to calculate the dielec- 
tric function spectra within the RPA employing norm- 
conserving pseudopotentials^. 



8tt 1 \ \ - p cv {k)p* cv (k) 

where indices c, v run over conduction (empty) and va- 
lence (filled) electron states, respectively. The spectral 
function F cv (ui,k) at the incident light frequency uj is 
given bjff^ 
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The normalization quantity £7 has dimension of vol- 
ume in direct space. The e represents the local part of 
the dielectric function. In this study we did not include 
non-localities due to the local field effect. Note that Eqs. 
([2]) and ([3]) represent an effective value for the dielectric 
function averaged over different tensor component of £ a ,fi 
(where a, f3 = x, y, z) in direct space. A more detailed 
explanation can be found iiPS 

The stability of various equilibrium configurations is 
judged by the adsorption energy. The adsorption energy 
per molecule is defined as the difference between the total 
energy of the adsorption system and the energies of the 
isolated components, namely the clean substrate and the 
adsorbate, divided by the number of adsorbed molecules 
N per unit cell: 



E a ds — 



Etot {E substrate ~1~ N 'E a dsorbate) 

N 



(4) 



A. Equilibrium geometries 

Table|l]lists the adsorption energies E a d s , per molecule, 
as calculated by Eqn. [4] and the optimized structural pa- 
rameters of the adsorbates on the surface. Atomic re- 
laxation of water on transition metals has be en studied 
both theoretically and experimentally^^ 22 1 23 1 . Two dif- 
ferent atomic geometries of water molecules physisorbed 
on Ag (111) surface are predicted: with hydrogen up and 
down. Our predicted geometries of physisorbed water 
molecules on Ag (111) surface agree, for the most part, 
well with those reported in Refs.l3"landl2"2"lwhere a similar 
theoretical framework was used. The H 2 0-up and H 2 0- 
down geometries are shown in Fig.[T^, and b, respectively. 

We also studied the effect of water molecule dissocia- 
tion on optical properties of silver. At room temperature 
the reaction H 2 ^ HO~ + H + occurs spontaneously. 
The proton is then free to enter into reaction with the 
surface. For a single water molecule adsorbed on a closed- 
packed Ag (111) surface the distance between the oxygen 
and the host Ag atoms, the O— Ag distance, of 2.78 A was 
reported in Ref. |22l Once enough H 2 molecules are 
added to the system to form hexamers on the surface, 
the molecules form an ice-like structure. For a single bi- 
layer of H 2 adsorbed on the surface, the zig-zag config- 
uration of ice isn't as pronounced, as for water molecules 
not in contact with the surface. Our oxygen-silver surface 
distances can be found in Table |TJ The water hexagons 
are slightly irregular as can be seen in Figs. [l^,b. The 
shift is much more pronounced in the H 2 0-up configura- 
tion. This asymmetry of the hexagonal structure of water 
molecules was previously reported for Ru{0001p3. Once 
one of the hydrogen is desorbed, see Fig. [TJ;, the remain- 
ing OH molecule shifts much closer to the surface and the 
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FIG. 1. Top and side views of the optimized unit cell configurations of water (a, b, and c), methanol (d and e), and ethanol 
(f and g) molecules adsorbed on Ag (111) surface. Panels represent the following configurations: (a) H 2 0-up, (b) H a O-down, 
(c) H 2 0-down-dsc, (d) CH 3 OH, (e) CH 3 OH-dsc, (f) CH 3 CH 3 OH, and (g) CH 3 CH 2 OH-dsc. The dissociate (dsc) configurations 
all have one hydrogen atom desorbed from the respective molecule and enter into a reaction with the surface. 

TABLE I. Chemisorption and physisorption adsorption energies (as calculated by Eqn. |4| and optimized structural parameters 
of H 2 0, CH 3 OH, and CH 3 CH 2 OH adsorbed on the Ag (111) surface. The average surface height is used for the O— Ag and 
H— Ag distances, representing the distance between the surface and the oxygen/hydrogen. Since there are two H 2 molecules 
per unit cell, two distances are provided. AAg is the difference in the average height of the top silver layer compared to a bare 
surface. O— H, O— O, and O— C represent the interatomic distances in the system; multiple distances are provided as needed. 



Structure 


E ads (eV) 


O-Ag (A) 


H— Ag (A) 


AAg (A) 


O-H (A) 


o-o (A) 


O-C (A) 


H 2 0-up 


-0.9047 


2.69,2.95 




+0.0007 


0.978,0.996,1.01 


2.78,2.95,2.98 




H 2 0-down 


-0.9471 


2.64,3.18 


2.06 


+0.0155 


0.994,0.998,1.01 


2.84,2.92,2.94 




H 2 0-down-dsc 


-0.1876 


1.83,2.25 


0.846 


+0.0233 


0.986,0.991,1.06 


2.51,2.94,3.37 




CH 3 OH 


-0.6017 


2.51 


2.24 


+0.0017 


0.985 


5.00 


1.42 


CH 3 OH-dsc 


+0.0622 


1.46 


0.968 


+0.0375 


2.94 


5.00 


1.39 


CH 3 CH 2 OH 


-0.6694 


2.61 


2.34 


-0.0026 


0.984 


5.00 


1.42 


CH 3 CH 2 OH-dsc 


+0.0004 


1.48 


0.911 


+0.0272 


2.96 


5.00 


1.40 



water hexagonal structure becomes even more irregular. 

Though the average surface height doesn't change very 
much, see AAg in Table |TJ the individual silver atoms 
on the surface do move quite a bit more. This can be 
seen clearly in the side view of Fig. [TJj. Another ob- 
servation is that the water molecules that do not lose 
their hydrogens are not oriented parallel to the surface, 
they are tilted with an angle of 7.69°. We report a neg- 
ative adsorption energy of -0.1876eV, as predicted by 
DMol 3 , indicating the stable structure. This is contrary 
to what was reported in Refs. and Previously, 
we carefully compared equilibrium geometries pred icted 
by CASTEP and DMol 3 with different moleculesPlSil 
and showed that the use of localized basis in DMol 3 



is beneficial in terms of prediction of stable geometries. 
Consequently equilibrium molecule-surface distances pre- 
dicted by LDA DMol 3 method better agree with those 
generated by VASP package including van der Waals 
contributions^ and agree much better with experiment. 

Methanol (CH 3 OH, Fig. [TJi) and ethanol 
(CH 3 CH 3 OH, Fig. fit) adsorb very similarly on the 
surface of silver (see Table III). Once the O-H bond is 
broken, the molecules rotate to bring the oxygen closer 
to the surface, however, the molecules also shift on the 
surface and position themselves such that the oxygen is 
between three Ag surface atoms. It is important to note 
that although a stable configuration was found with the 
O-H bond broken for both methanol and ethanol, those 
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configurations are metastable at best. Both methanol 
and ethanol (CH 3 OH-dsc and CH 3 CH 2 OH-dsc) have 
positive adsorption energies. One can speculate however, 
that the van der Waals force correction may remarkably 
modify the adsorption energies^ thus yielding stable 
configurations for dissociated configurations of both 
methanol and ethanol. 



B. Dielectric function and optical losses 

Optical losses of a system are proportional to the 
imagin ary pa rt of dielectric susceptibility function, 
(im(e)) ^ 10 l 26 [ Here we discuss the Im(s) spectra cal- 
culated according to Eqs. [2] and [3] for 7 monolayer 
thick Ag (111) slab with and without different adsorbed 
molecules. We also compare these results with the ex- 
perimental data for bulk AgpH. 

First, the predicted optical spectra for bare Ag (111) 
surface are considered. Without quasi-particle correction 
the predicted Im(e) LDA data calculated for a slab with 
experimental lattice constant shows a remarkable red 
shift in comparison with experiment. The LDA method 
used in this work produces compressed lattice of bulk Ag. 
Effect of lattice compression on electron energy structure 
results in blue shift of optical spectra^ thus improving a 
comparison with experiment as demonstrated in Fig. [2] 

Optical excitations from the surface states are respon- 
sible for the calculated increase of Jrn(e) of Ag-slab in 
visible and near infrared optical range with respect to 
bulk (see Im(e) data between 1 to 2.5 eV in Fig. [2]). 

In order to understand the nature of predicted changes 
in the spectra of Im(e) the projected density of states 
(PDOS) spectra are calculated for the systems studied. 
Details of the calculations of PDOS are described in 
Ref. Uni Our previous analysis of calculated band struc- 
ture and PDOS of the Ag-slab (without molecules in- 
dicated substantial contributions of surface electron en- 
ergy states which are in resonance with the bulk electron 
states and caused redistribution of the resulting density 
of states. In this paper we present the PDOS spectra 
analysis and focus on the effects on PDOS caused by the 
molecular adsorption. 

Optical transitions between the occupied states and 
Fermi surface are responsible for the predicted substan- 
tial increase of the Im(e) spectra of Ag (111) surface as 
compared with bulk in the near infrared region as demon- 
strated in Fig. [2] The transitions between the occupied 
states near —1.6 eV and Fermi level are responsible for 
the calculated shoulder at 1.55 eV in Im(e) spectra of 
bare surface, as shown in Fig. [2] 

Crystal truncation results in a predicted redistribu- 
tion of charge between Ad- and 5s-orbitals causing effec- 
tive energy shift of the 4d-orbitals towards lower energies 
(see Fig. [2]). Based on the PDOS data analysis we can 
state that predicted modifications of the Im(e) spectra 
of Ag (111) surface with respect to bulk are caused by 
surface-to-bulk and bulk-to-surface optical transitions. 
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FIG. 2. (a) Projected density of states of bulk (solid blue) and 
7 monolayer thick slab (dotted black) of silver. The shift in 
the d-orbital is due to the redistribution of electrons, (b) Cal- 
culated spectra of immaginary part of dielectric susceptibil- 
ity function for bulk (solid orange) and 7 monolayer thick 
Ag (111) slab (dotted blue) in comparison with experimental 
data measured in Ref. [27] (circles). 



This conclusion is in the alley of previously reported 
analysis^!, thus confirming bulk-to-surface nature of rel- 
evant electron transitions. Optical features predicted on 
Ag (111) above 3.5 eV are characterized by bulk elec- 
tron transitions with minor contributions of the surfaces 
electrons. 

Molecular adsorption substantially modifies spectra of 
Im(e). In this work we study the effect of molecular dis- 
sociation on the Ag (111) surface. At room temperature 
most solutions contain dissociated molecules which could 
interact with the surface states^. In Fig.[3]we present the 
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Energy (eV) Energy (eV) Energy (eV) 

FIG. 3. Projected density of states spectra of H a O (a), CH 3 OH (b), and CH 3 CH 2 OH(c) adsorbed on Ag (111) film. The 
thin (blue and red) lines depict the physisorbed configurations of (a)H a O-up (blue) and H 2 0-down (red), (b)CH 3 OH, and 
(c)CH 3 CH 2 OH; the dotted black line represents the chemisorbed configurations. The upper panel shows the total DOS and 
the other panel present orbital projected DOS as indicated by the panel indexes. The PDOS of bulk and silver surface can be 
found in Fig. [2] for comparison. 



effects of molecular adsorption on the calculated PDOS 
spectra of a Ag (111) film. On the bare Ag (111) surface 
the electron energy structure is substantially modified as 
compared to the bulk Ag. Molecular adsorption on the 
crystal surface and atomic structure relaxation resulted 
in substantial overal increase of DOS related to the 5s- 
and 4p-orbitals as shows in Fig. [3j We observed at least 
two pronounced PDOS 5s-peaks located near —0.4 and 
— 1.6 eV (occupied states) on Ag (111) surface as well 
as energy increase of both 5s- and 4p-unoccupied states 
in the region between and 4 eV. No substantial dif- 
ferences between surface and bulk of the 4d-unoccupicd 
PDOS has been observed. 

Analysis of the PDOS spectra shows that changes of 
molecular orientations on Ag (111) surface result in mod- 
ifications of mostly unoccupied p-electron states around 
1.6 eV and higher. The Im(e) function spectra calculated 
in visible and near infrared spectral regions are shown in 
Figg] 

In order to understand the effect of molecular dissoci- 
ation on optical losses of the Ag (111) surface we calcu- 
lated molecular configurations after the dissociation of a 
hydrogen atom from all molecules considered here. Ge- 
ometry optimization study indicates that the dissociated 
H-atom tends to interact much stronger with the surface. 
The same effect for water adsorbed on Ag (111), was re- 
ported earlier in Refs. [3] and [22] According to our to- 
tal energy minimization study, the dissociated hydrogen 
atom moved to the interstitial region. The adsorbed wa- 
ter pattern remarkably shifts while retaining the ice-like 
hexagonal structure. Dissociation caused enhancement 
of the interactions between unpaired oxygen and hydro- 



gen electrons. In all cases the dissociation of hydrogen 
leads to further substantial increase of Im(e) function 
which means the increase of optical losses of Ag (111) in 
relevant spectral regions. 

Comparative analysis of the calculated PDOS spec- 
tra of Ag (111) surface with dissociated molecules and 
that of bulk Ag show strong modifications of the elec- 
tron density: appearance of new occupied states located 
at around —1.7 to —2.5 eV below Fermi energy and un- 
occupied states near 0.5 to 0.7 eV above the Fermi level. 
Consequently in dissociative configuration there are con- 
tributions to the Im(e) function of electron transition - 
between occupied surface states and Fermi surface and 
between occupied and unoccupied surface states. The 
corresponding energy separations are between 1.5 to 3.2 
eV. Contributions of the transitions between Fermi sur- 
face and unoccupied host states to the observed changes 
seem to be minor. Our results therefore indicate that 
chemical reactions between dissociated H- and O-atoms 
occur mostly with unpaired 5s- and 4p-surface electrons 
of the Ag-atoms. 

Predicted strong increase of optical losses in infrared 
and visible range apparently dominates over that calcu- 
lated in the near UV region (see Fig. [4J . According to 
the PDOS spectra the (i-electrons are involved in optical 
excitations in the region above 3.2 eV. In contrast to the 
infrared and visible region, the molecular adsorption in 
ultraviolet somewhat even reduces optical losses in the re- 
gion above 4.0 eV (see Fig. [4} due to the discussed charge 
redistribution in the near-the-surface region. These find- 
ings should be interesting for the optical engineering of 
the UV devices. 
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FIG. 4. Imaginary part of calculated dielectric function of wa- 
ter (a), methanol (b), and ethanol (c) adsorbed on Ag (111)- 
%/3 x y/3 surface. The thin, dotted and dashed lines represent 
the undissociated configuration, while the bold solid lines cor- 
respond to dissociated molecules on the surface, respectively. 



Analysis of the data presented in Figs. [2] and [3] in- 
dicates several features related to the effect of molecu- 
lar adsorption on the Ag (111) surface. Physisorption 



only slightly modifies the Im(e) spectra mostly due to 
minor changes of atomic geometries. The most sub- 
stantial observed modifications of the Im(e) spectra are 
caused by chemisorption: molecular dissociation and cre- 
ation of new inter- atomic bonds. Adsorption of molecules 
substantially enhances optical losses in infrared region. 
Comparison with PDOS spectra indicate redistribution 
of the charge associated with 5s and 4p silver orbitals 
(see Fig. [3}. 

Results of this study clearly demonstrate strong con- 
tribution of the surface states to the optical losses on 
Ag (111) surface. Adsorption of organic molecules 
(ethanol and methanol) result in further increase of the 
losses on Ag (111) surface. Nature of the observed in- 
crease of the Im(e) function of Ag (111) surface induced 
by molecular adsorption is substantial modifications of 
surface electron energy structure and dramatic enhance- 
ment of electron transitions between surface and bulk 
states. 



IV. CONCLUSIONS 

Equilibrium atomic structure, projected density of 
states, and imaginary part of optical dielectric function 
spectra of Ag (111) surface covered with water, methanol, 
and ethanol molecules were studied using first principle 
density functional theory. Strong contribution of the sur- 
face states on optical losses is observed. Substantial mod- 
ifications of optical functions of metallic nano-slabs in 
near infrared and visible spectral regions, caused by sur- 
face states and molecular adsorption is predicted. Our re- 
sults indicate increase of the imaginary part of dielectric 
function of the Ag (111) surface due to charge redistri- 
bution of the surface Ag-atoms accompanied by the en- 
hancement of electron transitions probabilities between 
surface and bulk states. 
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